load("data_flies")

library(lattice)

flies.df$province <- as.factor(flies.df$province)

xyplot(DIC_n ~ year | province, 
       group = province, data = flies.df,
#       index.cond = list(c(31:1)),
       layout = c(4,8), cex = .5,
       type = c("l"), lwd = 2, 
       scales = list(relation = "free", x = list(cex = .7), y = list(cex =.5)), 
       par.settings=simpleTheme(col = "black", cex = .5),
       strip = strip.custom(
         bg = "lightgrey", 
         par.strip.text = list(col = "black", cex = .7, font = 3)),
       main = "Distribution of DIC cases (``flies'') by province", 
       xlab = "Year", ylab = "DIC cases")